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Aifcimci: 


In  this  paper  we  analyze  the  correlation  coefficient 
between  members  of  ary  two  colonies^  out  of  an  infinite  ensemble 
of  colonies.  It  is  assumed  that  each  colony  has  the  same  number 
of  members,  that  migration  takes  place  between  the  different 
colonies  ,  and  that  there  is  a  constant  rate  of  mutation  in  each 
colony.  An  explicit  formula  is  derived  for  the  correlation  func¬ 
tion  and  the  long  distance  form  of  this  .function  is  derived. 

It  is  shewn  that  under  rather  weak  restrictions  on  the  pat¬ 
ter?.  of  migration  the  asymptotic  form  of  the  correlation  func¬ 
tion  is  characteristic,  of  the  dimension  of  the  model. 


1.  Introduction 


The  theory  of  evolution  postulates  the  development  of  a  species 
by  the  fixation  of  mutant  genes.  Since  the  doctrine  of  survival  of  the 
fittest  is  no  longer  recognised  to  hold  in  its  strictest  form  one  can 
understand  differences  In  a  given  species  in  terms  of  the  development 
and  viability  of  several  alternate  genes  for  the  same  characteristic. 
Clearly,  geographic  separation  between  two  colonies  of  the  saw  species 
would  tend  to  promote  genetic  differences  between  members  of  the  two 
colonies.  One  can  interpret  the  formation  of  separate  races  in  this 
light.  It  is  the  purpose  of  this  paper  to  analyse  a  mathematical  mod¬ 
el  which  supplies  quantitative  information  on  the  effects  of  distance 
and  dimensionality  on  genetic  differences  between  groups  isolated  by 
distance.  The  model  was  first  proposed  by  one  of  us  (M.  K. ),  /I_7,and 
a  discussion  of  the  biological  implications  of  this  work  will  appear 
in  Genetics. /2  7.  The  present  paper  is  devoted  strictly  to  the  mathe¬ 
matical  development. 

Wright  was  the  first  one  to  discuss  quantitatively  the  pheno¬ 
menon  of  isolation  by  distance.  He  considered  a  large  population  in 
which  there  are  several  groups,  and  analysed  the  inbreeding  coeffici¬ 
ent  of  a  subgroup  relative  to  the  large  population  by  means  of  his 
method  of  path  coefficients  ,Z3-J7 •  Malecot  has  calculated  correlation 
coefficients  for  a  continuous  model  roughly  similar  to  ours,  ' 
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Uhe  present  model  considers  an  infinite  number  of  colonies, 
each  colony  containing  N  individuals.  It  will  be  assumed  that  breed¬ 
ing  proceeds  by  generations,  and  that  migration  may  take  place  between 
colonies  in  a  msnner  to  be  specified  more  precisely  below.  We  will  al¬ 
so  allow  autatlons  to  take  place  at  a  finite  rate.  In  this  way  the 
composition  of  each  colony  will  change  randomly  due  to  these  three 
causes.  We  shall  calculate  the  correlation  coefficients  S(p(l)  p(i+J]|) 
where  p(l)  is  the  frequency  of  one  of  the  alleles  in  colony  1.  We 
Ignore  effects  due  to  the  finlteness  of  the  number  of  colonies,  and 
we  shall  only  be  Interested  in  the  equilibrium  state  of  the  system, 

1.  e.,  the  question  of  the  time  of  spread  of  a  gene  through  a  popu¬ 
lation  will  not  be  treated  here. 

We  shall  begin  by  deriving  the  difference  equations  which 
describe  the  stepping  stone  model  In  any  number  of  dimensions .  Then 
we  &oall  specialise  to  the  one,  two,  and  three  dimensional  cases  when 
migration  Is  allowed  between  adjacent  colonies  only,  finally,  we  shall 
discuss  the  modifications  which  must  be  made  when  longer  range  migration 
is  allowed.  It  will  be  shown  that  under  certain  conditions  restricting 
the  fen*4*1**  of  long  range  migration  the  fora  of  the  correlation  func¬ 
tion  at  long  distances  is  that  predicted  by  allowing  the  difference 
equations  to  be  approximated  by  a  differential  equation.  She  present 
work,  in  its  mathematical  details,  is  the  discrete  version  of  the  theory 
of  multidimensional  random  processes  recently  discussed  by  Whittle,  A7- 
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2.  jformflgtlon,  of  the  Modal 

Let  us  consider  an  array  of  colonies j  for  simplicity  we  can 
think  of  them  as  being  found  on  tbs  coordinate  points  of  a  cartesian 
grid  which  trill  be  extended  to  Infinity  In  all  directions*  She  coor¬ 
dinates  of  a  point  will  be  denoted  by  4  •  where  on¬ 

ly  the  cases  n  *>  1,  2,  or  3  are  of  say  interest,  these  corresponding 
to  a  linear  habitat,  a  two  dimensional  dwelling  area,  and  a  three  di¬ 
mensional  dwelling  area.  She  latter  is  proposed  as  a  nodal  tor  micro¬ 
organisms  In  the  ocean.  By  convention  we  assume  that  the  caeponents, 

,  range  over  all  integers,  positive  and  negative. 

She  gens  frequency  of  allele  A ,  at  site  4k  will  be  denoted  by 
t*t>  *  In  our  model  the  gene  frequency  at  any  site  changes  la  any 
generation  by  three  mechanisms.  these  are: 

1.  Rxchwnge  between  the  given  subgroup  and  any  other  subgroup. 

2.  Mutation,  or  equivalently  exchange  between  the  given  subgroup 
and  one  located  at  Infinity. 

3*  Random  sampling  of  gametes  in  the  process  of  reproduction. 

It  will  be  assumed  that  migration  end  reproduction  occur 
periodically,  allowing  us  to  deal  with  epochs  rather  than  with  a  con¬ 
tinuous  time  variable. 

For  simplicity  we  shall  set  up  the  equations  for  the  £  IfM 
In  the  one  dimensional  case.  Ihen  we  shall  Indicate  the  extension  to 
the  n  dimensional  case.  In  one  dimension  we  let  vnj  denote  the 
probability  of  exchange  between  the  colony  at  <£•  o  end  the  colonies 
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at  and- j  .  Since  we  shall  only  deal  with  a  synaetrlc  situation  the 
probability  of  exchange  between  «fe»  o  and  -ft.j  will  be  .  If  "f* <■*■) 

denotes  the  gene  fraction  of  A,  at  jr  at  a  given  epoch  and  '(*)  de¬ 

notes  the  gene  fraction  of  A(  at  ■*  at  the  next  epoch,  we  can  write 

1*'(*>*  (?(*+>)+  +1*- 1 J) 

*  ( -f  («♦»)♦- f  («-»))♦...  ♦  *n„f  f 

where  !f  (*)  Is  a  random  variable  which  represents  the  change  due  to 
randan  luu«r^"fl  of  guetes,  and  jT  la  the  expected  fraction  of  gene  A,  In 
any  of  the  colonies.  The  term  "uf  represents  the  effects  of  antatlm 
In  rhiunging  +(*)  .  If  it  is  assuaed  that  there  is  no  selection  of  gene 
A,  then  has  the  properties  with  respect  to  expectation: 


eu'o). 

awfc 


(2.2) 


where  1b  the  effective  (constant)  size  of  each  subgroup.  This  para¬ 
meter  is  defined  in  terns  of  Hm  the  number  of  sales  and  Hf/  the  nun- 
ber  of  feaales  by 

±  *  1  I  J_  ±  \  (2*3) 

*•  1  W  ^ ) ' 

Equation  (2.1)  can  fee  cast  into  a  simpler  form  by  defining  a  new  vari¬ 


able 

j>  L*)  •  -'jPC*)  -  f 


(2.V) 
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which  eliminates  the  tern  •*\mf  •  Furthermore  we  will  only  be  concern¬ 
ed.  with  the  equilibrium  case  so  that  we  may  put 

♦  -  *(*).  (2.5) 

Finally,  in  order  to  rewrite  E^.  ix.  t)  in  a  simpler  form  we  shall 
Introduce  a  shift  operator  S.  This  operator  is  defined  by  the  proper¬ 
ties: 

s-fio  •  F(i-i) 

Making  the  changes  indicated,  we  can  write  the  equilibrium  equation  for 


«?<*>  as 

40MO  -  [(•- S"”*)  4. 


3  ?(.«>)  + 

±(«) 

(2.T) 

or 

(*)  *  L-jot*) 

(2.8) 

where  L  is  the  operator 

L  *  (l-  Z  - 

(2.9) 

In  higher  dimensions  the  operator  L.  takes  on  a  more  complicated  look¬ 
ing  form  because  there  are  more  indices.  For  example  in  two  diaan- 
sior  >  with  a  rectangular  coordinate  systen  it  is  necessary  to  intro¬ 
duce  two  shift  operators  8^  and  Sg  haring  the  properties 

l(L+>t  j) 

Sv  f  j  )  *  f  ( w, 


(2.10) 
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The  general  form  of  L  Is  then 


-  r  *■ 

■*/}»• 


‘j 


*■  K*  *."*)  (if-  V;J  (2.U) 


where  the  prise  denotes  oadssipn  of  the  term  l  *  Sm  e  . 

In  this  paper,  our  principal  concern  will  he  the  calculation 


of  the  correlation  function  which  we  denote  hy  f(*)  and  which  will 
he  defined  hy  the  relation 


r»).  E 


pit) 

puy 


(2.12) 


where  ^(jj  is  the  unnoreallsed  correlation  function. 

In  order  not  to  be  burdened  by  an  excess  of  notation  let  us  return  to 
the  one  dimensional  case  to  derive  the  equation  satisfied  by  r(* J 
If  we  multiply  £^.(2.1)  hy  $10)  and  take  expectations  we 
find  for  4?a0  the  result 

,<«)•  (2.13) 

However,  it  is  shown  in  Appendix  A  that 

E  [l?C«)  LfW].  e  [l?f  («)f(0)]  ■  L*f(»).  {tM> 

Hence  r(tt)  satisfies  the  equation 

(l-L')  r(H)  »  **0.  (2*15) 

Otis  result  is  also  valid  in  higher  dimensions  since  the  argument  which 
led  to  this  equation  does  not  refer  to  the  xuatioer  of  dimensions.  In  the 
case  of  a(o)  we  must  take  into  account  the  expression  for  the  variance 
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given  in  .  The  value  of  ^(o)  is  therefore  to  he  toil *4  from 

r  —  i 

■at 


m  fO-f)- 
2NC 


(I-  LX)f<0)  • 

— ■ »e  om 

(2.16) 

3.  Exact  Solution  to  the  Equations 

We  now  turn  to  a  discussion  of  the  solution  to  Eqp.  (2.l4) 
and  (2.16).  A  necessary  and  sufficient  condition  that  a  set  of  mot¬ 
hers  [r„)  to  he  a  correlation  function  is  that  it  admits  of  a  spectral 
representation,  /~ 9  7 .  Hence  it  is  natural  to  assume  a  solution  of  the 
form 

vr 

rt-i)  *  (jL/  j- J 

and  try  to  determine  the  function  Ffds^  h„)  •  We  will  carry  out 
the  detailed  calculation  for  *v*|  and  indicate  the  generalisation  to 
higher  dimensions.  Substituting  Eq,.  (3-1)  with  /  into  Eq.  (2.13) 
we  find  that 


r 

|  F(*)  («-^)  tr**®4*  •  o?  o. 


(3-2) 


It  will  now  he  demonstrated  that 

('  -  0-)  u*^9  •  C*  -  n’tuig))  vr*-H9  (3.3) 

where  H  (***•)  is  an  ordinary  function  of  Cn|  .  For  simplicity 


let  us  set  *•%  »  •-  £  -  *>• 

i*' 


.  Then 


.  i  ^  l*4'-  S’M  £  J£<y*s*fc) 


m  iZ"  £  l*1’’  ♦ *  sr***  *  s”’’1 )  «***• 


(3.4) 


»•*  4«» 
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by  straightforward.  Multiplication.  But  the  result  of  operating  on 
toft  by  the  operator  in  parentheses  is 

(s4+i  ♦  sUi  *  v  s"L‘*  )  tr.  *e 

(3.5) 

ss  Ua  (M+i*  tm  Cut  *».-  j) 6  «-*»  *■  u*(W-i.-i)® 


»  Itn-ftO  j)8*  t** 


*  M  te-a  <t 9  cr>  19  w»  j® 

Thus  we  find  for  £  u»  4R  8  "the  expression 

•  «  o* 

L  u*  4)9  *  ^  ZL  (.>*;%■  u*«.®  *«#• 

kM  jio  4 

‘  (l  tp  j8  if"  •  H*(u»®  )  t** 


(3.6) 


uhere  »  2.  +*  i  «**  1  %  * 

j*» 

Replacing  the  operator  O'  by  Hv it«  6)  in  Eq.  (3.2)  we  find 

(  Pl*)0-  H\we  »))  M  Nt  d®  *  O,  4a  a  o.(3*7) 
This  is  clearly  satisfied  if  we  choose  F(ft-)  to  be 


F(®) 


»-*\w»0)  (3.8) 

The  general  solution  for  rfe)  can  therefore  be  given  in  the  form 

*»  4® 


r(a> 


k  [  «-« « 
■r  *  “ 


(«-•) 


-he  o 


(3.9) 


The  constant  C  can  be  fixed  by  the  requirement  that  r(0)  *  I 


This  procedure  yields  the  result  valid  for  all  *  % 

»  (”<*»«»*»  ,  ("LsahT-i-  *  - _ 1 

r(<)|  -  *r)a  1 1 - H(m»9)  ^ 

tTS 


*  *•  »  -  H\u»d) 


4  J _  “ 

HUM)  \  *  HCu*») 


10) 

4® 
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In  higher  dimensions  the  same  reasoning  leads  to  the  representation: 


j_  (  (  Cm  4ft.il..  A 

<“>'  J  I-  r3-U> 


tuy*  3  J 


4et  4»v-  *©« 


Finally,  we  return  to  the  calculation  of  p(0)  fron  Bq..  (2.16)  now 
knowing  the  values  of  the  r(*)  .  The  general  procedurt  for  finding 
this  quantity  is  to  substitute  p o)  K*)  for  the  value  of  p(*)  . 
In  this  way  we  find 

v  (3*12) 

p(o) 


which  is  also  true  In  higher  dimensions. 

The  value  of  0-  )r(0)  ±B  found,  from  the  representation  of  Eq. 

(3.11)  using  the  Identity  of  (3.6).  We  have,  in  fact 

0-  I?)  r(0)  «  L  . X *  i9  (3#13) 

5r  I  —  — 

•  I-  ^(cr*©) 

which  Implies  the  expression  for  f>(°)  • 

p(0).  tO-f)  ,  f  f)  <3-»> 

utM.o-L'ino)  “  4Htr 

- 

«  t  -  HxCw»  ») 


10 


When  <>*»  »  0  but  -wt  is  greater  than  zero  then 

I  V  " 


p(0)  - 


fO-f) _ 

,+  U-  **«,) 


(3.15) 


which  agrees  with  a  result  given  earlier  tor  Wright,  /%  7  •  She  result 
analogous  to  Eq,.  (3.10)  holds  also  in  higher  dimensions. 


4.  Specific  Solutions  to  the  General  Equations 

In  this  section  we  shall  present  a  study  of  several  specific 


examples  using  the  general  techniques  so  far  developed.  In  addition 
we  shall  derive  an  asymptotic  expression  for  r(£)  under  the  assump¬ 
tion  that  ♦  49*  is  large.  She  simplest  trample  is  that 

of  a  habitat  in  which  migration  can  accur  between  adjacent 

colonies  only.  Tor  this  case  only  %  differs  from  zero  and  the  func¬ 
tion  H(to  •)  is  given  by 


H(to»).  •*,  0- «•»•). 


(4.1) 


She  relevant  integrals  can  all  be  derived  from  the  single  formula 


1  f 

n  t: 


i.9(l0  1 


use 


-»)\ 


=  1-0 


I 


O'. 


With  these  results  we  find  for  r(<fc) 


x  >  i 


x  *  -  • 


(4.2) 


(4.3) 


r(*h> 


A,(«)  *  AXC«) 
A,U»  *  A t(o) 
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where 


0*  *  -  jMsjHsSr]*  (k,k) 

-  ’-^T 


Althou^i  the  solutions  so  for  have  been  perfectly  general,  there  are 


approximations  possible  in  writing  the  solution  based  on  the  fact  that 

<♦1^,  «  (4.5) 

for  all  realistic  situations.  In  fact  we  shall  generally  assume  that 
is  of  the  order  of  /©"'  and  'v*»i»  is  roughly  4*  ;0~ r  . 

This  implies  that  A,  (a)  is  small  in  comparison  to  A,(*J  .  This  is 
also  true  in  higher  dimensions.  For  example,  in  the  present  case, 
with  the  parameters  that  we  have  just  mentioned 


.  *  (4.6) 

A,(«)*  177  (  972) 

\  (H)  «  0. 2C1  1-M4)*  . 


It  will  be  noted  that  with  the  present,  general  method  of  solution, 
one  does  not  have  to  choose  the  solutions  with  the  property  r  l«)  -*  o 
This  comfit  ion  is  insured  by  the  trigonometric  representation  of  r(*) 
given  in  Eq,.  (3.11)  and  the  Riemann-Lebesgue  lemaa,  /~10_7 . 
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Let  us  next  consider  a  two  dimensional,  rectangular  lattice 
of  colonies  in  which  migration  can  take  place  only  to  the  nearest 
neighboring  colonies,  as  pictured,  in  Tig.  1.  So  be  completely 
general  let  us  suppose  that  migration  in  the  >  direction  Is  differ¬ 
ent  from  that  in  the  g  direction.  We  may  use  the  general  form  of 
Eq.  (3.11)  and  Introduce  the  function  of  two  variables  tf  («*•  8t>  %) 

ty 

a  I-*.-  m,(i- 

She  two  Integrals  which  figure  In  the  evaluation  of 
are: 

w  «r 

A, (%‘S.)*  .-v  (  1 

*  •  ’M'-U't,)*  (h.Q) 


*itx 


•  •  A-eu-a*, (i-u. *)-  *.v(i -Cn^j 


These  may  be  transformed  Into  single  Integrals  which  preserve  the  sym¬ 
metry  inherent  in  them  by  use  of  the  identity 


±- 

ft 


c-  **  « 


(4.9) 


to  put  the  denoolnators  in  an  exponential  form.  In  this  way  we  find 
that  ft,  *%*».)  and  Av(«,  r\v)  can  be  expressed  alternatively  as 


,(%<%)•  £  dt  . 


(4.10) 


n'r%,r  „  #  tu.9,  JQ  dt, 

•  e 
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»  <^v  ^  «~C1  "** |  j  u*n,»,  im  %•, « 

•  o  o 

(-.,tw.»,  <-*,ti.«l  . 

Bw  stop  of  interchanging  orders  of  integration  can  toe  Justified,  in 
detail.  She  angular  integrations  can  toe  carried  out  toy  eeans  of  ttoe 
Bessel  function  representation,  > 

£  »9  eT*  '*'•  a  9  *  I*  t*"*)  (*.U) 

sliers  3»  Is  a  Bessel  function  of  laaglnaiy  argunsnt.  In  this  say 
sa  find 

h, (., ,«j  •  i-  (V  <“**"'*  -u*  (».m) 

o 

\KM'  J't- 

o 

It  is  dear  fron  this  representation  that  for  the  values  of  ,  *» , 
and  a>% considered  here  A,  (*.,«%)  is  always  greater  than  the  ab¬ 
solute  seine  of  av)  since  X*(»j  is  posit  its  for  no  and  the 

quantity  *»„  «•  w*(*  is  such  las#  than  1  "S  *  wv  . 


Ik 


The  diagonal  terns  A, ,  (">  n )  can  he  evaluated  exactly  In 
terns  of  Legendre  functions  of  the  second  kind  or  in  terns  of  complete 
elliptic  integrals.  It  is  known  that 

*  ^Br  9-i 


so  that  if  we  replace  t  by  lb  and  C.  hy  it  we  obtain 

Where  the  argunent  of  the  Legendre  function  is  negative  since  a>  l»4  c. . 
for  <k*o  we  use  the  integral  respresentatlon  of  the  Legendre  Auction 
to  transform  this  into  a  fora  more  convenient  for  calculation: 


<V"»*  7£r  ( 


it 


ta 

.  Wi  I 


(*.15) 


it  Uil 


■  *  Ws ) 


where 

K(V).  (  -jg»  (“-l6) 

'«  'll-  ** 

is  a  complete  elliptic  integral  of  the  first  kind.  Thus  the  nozaalis- 

m 

ing  constant  C  is  found  from 

r  •  ,  r  1  v  ( _ ! - -  V>.i7) 


V' 


K  (— r 

M*  (*- 


U*- (x'^)  - 

dw*,e\,  V  ’  v 
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When  /*»,  and  r**  are  smell  In  comparison  to  f  we 

may  derive  a  simpler  expression  for  C"1 .  Let  us  consider  the  first, 
elliptic  integral  in  Eq,.  (4.17)  and  set 


€  * 


<***,¥*v 


(4.18) 


i.  e.j  «  is  also  small  in  comparison  to  I  .  The  argument  ()♦*)"' 
is  therefore  close  to  \  and  we  may  then  approximate  the  elliptic 


integral  ty 


xCrh)-  ou-*~«).  <*•*» 

In  the  second  elliptic  integral  to  appear  in  Eq.  (4.17)  notice 
that  the  term 


U-  *^»)V 

IW,  1WV 


is  large  in  comparison  with  i  so  that  the  argument  is  small, 
can  therefore  use  the  approximation 

*<>)>■  ?t>‘  1*  Jo'*  0lf>l 

Hence  an  approximate  value  for  C*'  is 


(4.20) 
We 

(4.21) 


C* 


I 

IT 


(4.22) 


-!■  _ 

2  \J  (*-*’•)'■-  »U- 


(x  -  "*«)j 


1 6 


We  have  neglected  tense  of  the  order  of  /  Art*  end  *%  / 
since  these  are  0 (io'k)  for  all  prohleoe  of  interest. 

Exact  results  can  he  obtained  for  the  diagonal  elements  by 
recursion.  The  expression  for  is 


L«-i  u.*0-£j 


(*.83) 


We  have  already  related  Q.^(-*)  to  the  complete  elliptic  Integral 
X  (  j  .  A  similar  calculation  shows  that 


Qt  u ) . ;  £uT7T  6  (£, )  -  *  t&, )  ("•») 


where 


EC*)  1*  the  complete  elliptic  integral  of  the  second  kind 


Values  of  far  **  >  I 

lation,  /llj 


can  be  obtained  from  the  recurrence  re- 


c- » i . a  t  [i.  - 1 1  o..  k  <-#)  - o..t  i-»  )J .  (l,-a6) 


It  does  not  seen  possible  to  derive  values  of  f  (%'S)  da  general  for 
>*,♦  in  terns  of  tabulated  functions.  However,  in  the  syaartrlc 
case,  r%  ,  it  is  possible  to  derive  expressions  for  p(%i%) 
either  by  appealing  to  the  difference  equation  or  by  returning  to  the 
integral  expression  for  the  A'*  .  For  example  to  get  ft  (o,  i)  we 
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note  that  when 

A (0,1).  AC‘,0)*  1  [A<o,0*  A  C »> oj ]  .  ^,27) 

Hence 

#  (4.28) 

A,  (0,0  ♦-  VM-  i(2*£»)A,‘(4°)  ♦ 


fron  vhlch  &  C©>  i)  can  be  calculated.  In  theory  it  would  be  possi¬ 
ble  to  build  up  a  couplets  set  of  the  *»J  for  the  sysuetric 

case  by  recurrence.  However,  this  quickly  becoaes  tedious  sad  is  in¬ 
deed  quite  special  because  of  the  restriction  .  Hence  we 

turn  to  a  different  range,  naaely  large  values  of  *  *t.  • 

tor  large  ft,*  another  approach  is  possible.  Let  us 
return  to  the  expression  of  A,|%  (ft,,  ftv)  1“  terns  of  trlgouoavtric 
Integrals  (Eq.  (4.8^.  We  can  write  the  expression  for  A,(ft„*g  •• 


«■  w 


V«„  rs.)  *  ,Tv  j  |  *"'*'*. 

•  •  *%»♦  m, (»-  a**#,)  ♦ 


(4.29) 


by  syaatry. 

A  slailar  reduction  in  the  area  of  integration  of  AvCft(Jftfc) 
also  holds.  In  the  Halt  *n>l;  r\x-+  m  we  notice  by  asynptotlc  ar¬ 
guments  slailar  to  those  given  by  Duffln,  /”l2 J,  that  the  aaln 
contribution  to  the  integral  cones  froo  the  neighborhood  of  1*9,,  ®v)  —  ( 0,  C  j 
since  the  dencnlnator  is  a  nlnlaua  there.  Hence  the  integral  can  be 
approxiaated  by  expending  the  denoolnator  around  (O,  0  )*• 
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J_  |  I 

^  v.  **♦  Ei^s 

^  1. 


(4.30) 


Again,  in  the  Halt  '*.1>  *v  — ♦  m  the  liaite  of  integration  can  he 

fixed  as  (o,*  )  with  an  exponentially  snail  error.  *»">■  we  find 

“  *• 

|  *«o  *iOt  (4.31) 

•  •  A\»  *  *1 8^  + 

Cr»  l*«  S,  *.)  J  <**,  **». 

»"».  «•  •»*  1- 


In  a  siail&r  way  it  can  he  shown  that  an  appxoxlaate  expression  for 


*tK,  *v) 


is 


fj 


.  0*-3£) 

*("*,*  r\))  *7C^"t 


I 


Both  of  these  integral  depend,  aside  froa  the  trlgonoaetrio  factors, 
only  on  .  Bochner  has  shown,  [ l]J ,  that  if  an  Integral  is  of 

the  forn 


UK/*, i‘"J 
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depends  only  on  the  quantity  •<  *  C*N***'''  ♦o£)'/* 


and  the  integraJ  can  be  reexpressed  aa 


(4.3»0 


where  J«-y  Ct)  is  the  Bessel  function  of  order  # 

In  the  present  case  we  nay  set 


« .  ( £  ♦  *2.  \k 
v  w*.  w.  / 


and 


(*•35) 


(^.36) 


'"'V- 1 

by  /jj/vhere  K,U)  is  a  modified  Bessel  function  of  the  third  kind. 
In  a  similar  way  we  find  that  Aj(y»,}av)  is 

.  ■*,♦  *v _ .  (A.3T) 

4T\Pvwv 

The  parameters  and  are  of  the  order  of  magnitude  of  S*  i o“K 
and  r*m  is  roughly  k*\0“r  •  Hence  if  at,  and  7\  are  both  9  or 
greater  then  ti  is  greater  than  70  and  the  argument  in  the  Bessel 
function  approx  lost  ion  for  A,  *v)  is  roughly  6 
argument  in  the  expression  for  A^^rtJ  is  over  loo. 


while  the 
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Considering  that  the  function  K0  («)  has  the  asymptotic  expansion 


*e 


(4.38) 


we  see  that  Avt*ijA^  is  negligible  in  c  caparison  with  A,  (*,,wv) 
and.  that  for  one  can  satisfactorily  write  for 


•U*} 

m, 


2  *v*»  n*  \" 
—  )  ' 


(4.39) 


In  particular  if  ^|»wivand  we  set  >/n^v  »  R 


rl*) 


ifair 


,  we  find. 
(4.40) 


Thus  we  see  that  in  contrast  to  the  one  dineosioaal  case  in  ***< the 
correlation  falls  off  exponentially  at  large  distances,  the  two  diswn- 
sional  correlation  function  falls  off  as  an  exponential  divided  by  JW  . 


It  olgbt  he  objected  that  such  a  conclusion  depends  critically  on  our 
model,  1.  e.,  that  If  longer  range  migrations  are  included  the  asymp¬ 
totic  form  of  the  correlation  function  would  he  changed.  In  the  next 
section  we  shall  show  that  if  the  long  range  migrations  axe  sufficient¬ 
ly  weak  the  general  form  of  the  correlation  function  will  rams  in  un¬ 
changed  except  that  v*,  and  era  replaced  by  different  combina¬ 
tions  of  the  w\'t . 

A  closer  analysis  shows  that  the  next  tem  in  the  erpanston 
is  of  the  order  of 


ft 


Which  is  ' 1’lm*  than  the  result  of  Bq,.  (4.40)  by  a  factor  of  Jr 

Finally  we  may  turn  to  the  three  dimensional  ease  in  which 
migration  is  characterized  by  three  exchange  coefficients, 
and  Wj  as  shown  in  Fig.  2.  For  this  case  we  use  the  same  reason¬ 
ing  as  in  two  dimensions  and  thereby  find  the  results 


0  . 

A-KSMs  -i-r  (T( 

T  i  n,  v  a 


la  this  case,  an  expresnioa  ia  closed  form  for  p(ot0{0)  is  known  only 
for  »*»•  -  0  and  ** ,  *  Mv  *  "»j  ,  otherwise  numerical  nethods  must  be 
used.  Fortunately  tables  are’ aval labia  which  axe  useful  for  the  case 
**>  •  **«.  IP\J*  Ihis  is  the  case  of  Isotropic  migration  la  a  plane 
but  a  different  aaount  of  migration  la  the  third  distension.  Such  a 
■odel  would  be  useful  to  describe  the  composition  of  aa  oceanic  habi¬ 
tat  la  which  it  sight  be  supposed  that  there  is  no  distinction  be¬ 
tween  the  two  directions  of  horizontal  migration,  but  there  is  a  dis¬ 
tinction  between  migration  In  a  horizontal  or  a  -vertical  direction. 

The  tables  in  [}hj  axe  expressed  In  terms  of  the  function 


ItSV-S-  if)  •  fc  1  S I  *»'»<■» 


(V.te) 


•  ~  C*.  X  -  Or.  —  AC«> 

Performing  the  necessary  redactions,  we  can  express  the  A'*  In 
of  this  triple  Integral: 

<*.*3) 

At  V 

■  t-0  ^  !*,.*>-  i"i  t  w>>  .  1~nn*  -  >  ). 

1  *S  *  a«y**»,  ' 


The  -value  of  plo,o;o)  Must  be  obtained  froai  the  tables. 

^  is  such  greater  than  i  ,  the  -value  of 
is  given  by  the  asymptotic  expansion: 


When 
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pK.-NV-  |/|  4tR  «•' 


Aft 


(4.V0 


where 


R'.  v  **, 

A1-  a  U*  W  — V  \  . 

'  V  / 


(4.A5) 


In  particular  if  m  •  then  the  correlation  function  has  the  charac¬ 
teristic  fbrm 

Rte-  (*•*> 

4irR 


where  Rv.  y»J*  «.  *»  .  This  is  seen  to  flail  off  such  acre 

quickly  than  the  one  or  two  dimensional  cases.  Corrections  to  Eq,. 
(4.46)  axe  given  in 

In  Figure  3  we  have  plotted,  for  the  purpose  of  conparlaon 
some  results  for  one,  two,  and  three  distensions.  In  these  we  have  set 
r*m  *  4s  t  o"x  and  all  of  the  other  rw '«  are  set  equal  to  one  an¬ 
other.  In  one  dimension  we  have  put  *\»  o.i  ,  in  two  dimensions 
ry1) .  a. or  and  in  three  dimensions  <*»»  *  o .  o  7  I  2  . 

It  can  he  seen  that  an  increase  in  the  number  of  distensions  has  a  crit¬ 
ical  effect  on  r  <  )  . 


2k 


5*  general  Interactions 

So  far  we  have  restricted  ourselves  to  a  rectangular  lattice 

with  migration  permitted  between  nearest  colonies  only.  It  might  be 

thought  that  the  conclusions  that  we  have  drawn  concerning  the  ie«g 

range  form  of  the  correlation  function  depend  on  the  particular  modal 

we  have  chosen,  and  that  If  longer  range  migrations  are  allowed  our 

conclusions  are  Incorrect.  Ve  shall  show  that  if  the  long-range  mi- 
>» 

gration^ sufficiently  week.  In  a  sense  to  he  specified  more  exactly, 
the  correlation  function  retains  the  same  form  for  long  distances  as 
that  found  for  near  colony  migration.  For  simplicity  wa  ahall  derive 
the  result  in  one  dimension  and  indicate  the  extension  to  higher  di¬ 
mensional  habitats.  She  reasoning  for  these  rnmslnn  exactly  the 


In  one  dimension  the  representation  of  Eq,.  (3. 10)  Is  val¬ 
id.  For  large  *1  ws  expand  A,(M  and  A,(M)  as  Indicated  in  Bq. 
(4.30).  For  A,  l")  we  find 


....  i.  I  a*.  0  dS 

A,  O')  XT  j  - 

'0  *■  m,(i-  Lr»0,)  + 


—  2-  (  C«  d$ 
XT  1  * - J* 


'0  %  ♦  2.  T  i* 
a 

1*  ’ 


(5.1) 
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where  we  hare  defined  the  paraaeter  by 


j*' 


(5-2) 


Urns  it  can  he  seen  that  under  the  hypothesis  that  ^  i e  finite 
A,  decreases  exponentially  with  n  ,  the  coefficient  of  w. 

being  -  (2 r*+  “  waa  *oun4  earlier  for  the  special  ease 

of  nearest  colony  nigraticn.  She  function  A*(w)  for  large  n 
is  given,  in  the  sane  Banner,  by 


in  which 

m 

*  T  **i 

i* ' 


(5.4) 


It  is  necessary  to  find  the  value  of  p(o)  by  sane  nueerical  or  approx- 
iaate  method  in  the  present  instance,  notice  that  the  inequality 

I  >  (5.5) 

always  implies  that  A,(**) b  I  Av(«)  I  ,  or  that  the  first  term  is 
always  doninant  in  the  liait  r>  -»  m  .  However,  the  possibility  of 
long  range  migration  can  make  the  contribution  of  l*»)  signifi¬ 
cant  in  the  range  of  n  which  nay  be  of  interest. 
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In  two  dimensions  we  oust  consider  an  array  of  quantities 
in  which  is  the  percentage  of  the  colony  at  (o>Oj 
which  is  exchanged,  with  the  colony  at  (ly  j )  .In  terns  of  these 
quantities,  and  auxiliary  functions  defined  by 

m  m 

A,,'  z  Z  irJS**ij  (5.6) 

J  '  •  1*0 

we  find  In  two  diawns  iocs  that 

~ _  K,  (  a  \TT*m  )  (5.7) 

Aj.Kj'HW  (-  »  "v  K0  0*  J  3  U-iea-^a"  *Aes  J  ) 


where 

h‘^o 

Thus  we  see  that  the  only  difference  between  results  derived  for 
nearest  colony  Migration  and  for  the  present  ^restricted  case  is 
that  and  mv  are  to  be  replaced  by  f*.^  and  respectively. 
When  y-t0  and  are  not  finite  this  account  is  not  valid  and 
other  fores  are  possible  for  r  *%. )  • 


.  (5.6) 

I!*  . 

Na 


However  when  migration  can  take  place  between  a  finite  group  of  colon¬ 
ies  only,  the  asymptotic  form  of  the  correlation  function  will  always 
be  as  given  in  Eq.  (5.7)  *  She  conclusions  for  the  three  dimensional  case 
case  axe  also  the  same. 

Results  for  different  periodic  arraagaaents  of  colonies 

1 

resemble  those  given  in  this  paper,  formal  expressions  can  be  derived 

I 

for  r(£)  in  terms  of  mltiple  trigonometric  Integrals  as  la  the  pre¬ 
sent  analysis.  She  conclusions  which  can  be  drawn  concerning  the  long 
distance  form  of  axe  the  same  as  given  In  the  present  paper, 

and  aipie*  moat  elgnlflcaatly  on  the  il  Inins  inn!  itj  of  the  habitat, 
■atuarally,  the  valma  of  r(4&j  for  small  *  varies  between  different 
types  of  lattices  and  must  be  calculated  Individually.  One  can  write 
down  a  continuous  version  of  the  present  model.  She  analysis  leads 
to  partial  differential  equations  analogous  to  our  partial  difference 
equations,  and  Indeed  the  same  farm  of  ri* )  for  large  A  can  ha 
recovered*  However,  in  two  and  three  dimensions  the  solutions  to  the 
partial  differential  equations  are  singular  at  the  origin  and  ad  hoe 
methods  axe  required  to  fix  the  normalising  constant.  It  Is  therefore 
felt  that  the  present  model  is  a  more  accurate  representation  of  zeal* 
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Appendix  A:  Proof  of  the  Equality  E  [l*(«0  If  (0)]  *  E  W  $ 

In  this  appendix  we  prove  the  equality  stated  in  the  title 
for  the  one  dlawnsiomal  case.  A  staple  extension  of  the  notation  suf¬ 
fices  to  prove  the  result  for  higher  dimensions.  The  proof  is  straight¬ 
forward  and  consists  of  a  direct  calculation  of  the  terns  Involved. 

We  my  write  for  the  left  side  of  the  equation 

£  «  i6  [I  S"'r)f(0)] 

«  [l  I  **' r  (f <-0*  ?(«•))] 

'  u  jt0  rto  *• 

*  J  £  f  i  2  «*  f (-rif(afj) 

*  •-  jttf  r.a 

t  j  J  +  fl-r)  -p  l«-jj  ] 

m  - 

*  J  I  2  f  p(.«tjrr )  «-  P  (*-  l“r  ) 

“  j*o  r*o  1 
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For  the  right  side  we  find 

£ *  E  tf<°)  L\f (*)] 

-  i  t  Lf'°l  i  £  *-.»>,  lS*>Vj)(S'r 

M  j*a  N, 

n  oo 

r  5*  £  X  tf\r  (  p(*M*  r)  ♦  <*)  *> 

*  c  Ll  f  i«)  L£(o)]  . 

Another  proof 'Of  this  equality  can  he  derived  hy  a 

f<*).  EJ 


V')  f  <*)  ] 

lens  of  the  aoeotral 


representation  of  . 
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Figure  Captions: 

Fig.  1.  Plot  of  r«j  u  a  function  of  for  one,  two,  and 
three  linens Ions . 

Fig.  2.  Schematic  representation  of  exchange  between  colonies 
In  a  two  dlnanslonal  habitat . 

Fig.  3.  Schanatic  representation  of  exchange  between  colonies 
In  a  three  dimensional  habitat . 
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